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i[r We study the motion of an inertial particle in a fractional Gaussian ran- 

'~~; doni field. The motion of the particle is described by Newton's second law, 

1 -^ where the force is proportional to the difference between a background fluid 

Cd velocity and the particle velocity. The fluid velocity satisfies a linear stochas- 

^ tic partial differential equation driven by an infinite-dimensional fractional 

Brownian motion with arbitrary Hurst parameter H £ (0, 1). The usefulness 

y—i of such random velocity fields in simulations is that we can create random 

^' velocity fields with desired statistical properties, thus generating artificial 

.^ images of realistic turbulent flows. This model captures also the clustering 

OS phenomenon of preferential concentration, observed in real world and nu- 

^^ merical experiments, i.e. particles cluster in regions of low vorticity and high 

^£^ strain rate. We prove almost sure existence and uniqueness of particle paths 

fvj and give sufficient conditions to rewrite this system as a random dynamical 

1—1 system with a global random puUback attractor. Finally, we visualize the 

h>. random attractor through a numerical experiment. 
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1 Introduction 

Motion of inertial (i.e. small heavy) particles in turbulent fluids occurs in natural 
phenomena as well as in technological processes and therefore has excited theoretical 
investigations, see for example the book [46] and the references therein. Examples for 
such processes are formation of raindrops, evolution of clouds and combusting of liquid 
fuel. 

The starting point for many theoretical investigations concerning the motion of inertial 
particles in turbulent flows is Stokes' law (see e.g. [3]), which says that the force exerted 
by the fluid on the particle is proportional to the difference between the background 
fluid velocity and the particle velocity, i.e. we are concerned with the following transport 
equation 

Tx{t) =v{x{t),t) -X{t), (1.1) 

where v{t, x) is the velocity of the fluid at point x in space at time t and x{t) is the 
position of the particle at time t. Here the response time r = ^ in two dimensions or 
r = Q^ in three dimensions is often called Stokes ' time, where m is the particles mass, r 
the particles radius, u the fluid viscosity and C > a universal constant. We neglect that 
C actually depends on the radius r and the relative velocity, making the law non-linear 
[3]. Further, an important non-dimensional parameter related to the equation (1.1) is 
the so-called Stokes' number St, which is the ratio of the particle aerodynamic time 
constant to an appropriate turbulence time scale. In turbulent fluid flows St is usually 
defined by St = t/t^, where r^ = e~^''^u^''^ is the eddy turnover time associated to the 
Kolmogorov length scale rj = e~^/^i/3/'^ with viscosity v and mean energy dissipation rate 
e. If r — )• oo the equation (1.1) tends to x(f) = 0, which is the equation of motion a of 
particle moving with constant velocity. And if r = 0, i.e. the inertia of the particle is 
neglected, we get x{t) = v[x{t),t), which is the equation of motion of a fluid particle 
or a passive tracer, which is a particle that follows the streamlines of the fluid. Various 
extensions of the basic model (1.1) have been considered in the literature, in particular 
by Maxey and collaborators [24, 25, 45]. 

Real world and numerical experiments show that the distribution of inertial particles in 
a turbulent fluid is highly correlated with the turbulent motion ([13, 14, 36, 42, 43]). 
The particles cluster in regions of low vorticity and high strain rate. This clustering phe- 
nomenon is known as preferential concentration. The Stokes' number St plays a central 
role in the effect of preferential concentration. If the particle and fluid time constant 
have almost the same order, i.e. St ^ 1, the particles concentrate in regions where 
straining dominates vorticity. Experiments at high or low Stokes' numbers do not show 
this clustering phenomenon. 

A model for the motion of inertial particles in two dimensions which covers the pref- 
erential concentration phenomenon was introduced by Sigurgeirsson and Stuart in [40] 
and analyzed in a series of papers [16, 21, 30, 31, 32, 33, 34, 41]. This model consists 
of Stokes' law (1.1), where the velocity field is a Gaussian random field that is incom- 
pressible, homogeneous, isotropic and periodic in space, and stationary and Markovian 



in time. This gives the equations in non-dimensional form 

Tx{t) = v{x{t),t)-x{t), (x(0),i;(0)) gT^ X M^, (SI) 

v(x,t) = V^i;{x,t) = (^^^{x,t),-^^{x,t)), {S2) 

dipt = jAiPtdt + -f'2Q-2dWu i^o^^V, t > 0, {S3) 

where r, 7 > 0, T^ is the two-dimensional torus, W is an infinite-dimensional Brownian 
motion with an appropriate infinite-dimensional separable Hilbert space V as state space 
and the self-adjoint operators A and Q^onV will be chosen to match a desired energy 
spectrum of the velocity field. The equation (53) is interpreted as a linear stochastic 
evolution equation ([12]). Precise assumptions will be given later. In [40] various quali- 
tative properties of the system (S) have been studied, such as existence and uniqueness 
of solutions and existence of a random attractor. In particular, numerical simulation 
in [40], see also [6], indicates that system {S) also covers the preferential concentration 
phenomenon, where now the parameter r can be interpreted as Stokes' number, whereas 
7 indicates how fast the velocity field decorrelates. Therefore, the random attractor of 
the system (S) is highly relevant for the study of preferential concentration of inertial 
particles. 

Further, various limits of physical interest in system (S) have been studied: rapid 
decorrelation in time limits [21, 30, 31] and diffusive scaling limits (homogenization) 
[16, 32, 33, 34]. 

Also of interest are studies of system (S) for fluid particles or passive tracers, i.e. r = 
and so (5*1) is replaced by x{t) = v{x{t),t). This problem in a similar framework was 
considered among others by Carmona, see [7] and references therein. 
In this work we generalize the results of Sigurgeirsson and Stuart in [40] to the fractional 
noise case, i.e. in the system (S) the stochastic evolution equation (53) will be driven by 
an infinite-dimensional fractional Brownian motion B^ with arbitrary Hurst parameter 
H G (0, 1). We recover system (5) as a special case for H = ^. In particular, we improve 
Sigurgeirsson and Stuart's assertion in [40] concerning the existence and uniqueness of 
the random attractor in such a way that we extend the universe of attracting sets from 
deterministic bounded sets to random tempered sets with an explicit representation of 
a random tempered the universe absorbing set. The main motivation to use fractional 
noise, to model now a fractional Gaussian random velocity field, is that there exist cer- 
tain statistical similarities in the scaling behaviour of a fractional Brownian motion and 
a turbulent velocity field. These statistical similarities are described in the next section. 
Also of interest and what Sigurgeirsson and Stuart have not done, is to give analytic 
(upper) bounds of the Hausdorff dimension of the random attractor, which is almost sure 
constant due to the ergodicity of the noise. This will be studied also in the fractional 
noise case in a forthcoming paper. 

The remainder of this article is organized as follows: In section 2 we present the mo- 
tivation to use fractional noise to model the random velocity field. Section 3 provides 
properties of a stationary fractional Ornstein-Uhlenbeck process, which are needed in 
Section 4, where we introduce the generalized model (S) for motion of inertial parti- 



cles in a fractional Gaussian random velocity field and prove almost sure existence and 
uniqueness of particle paths. Section 5 describes how to match desired statistical prop- 
erties, in particular a favoured energy spectrum of the velocity field. Further, we verify 
in Section 6 that the system (5) defines a random dynamical system, which admits a 
random pullback attractor. Finally, we visualize the random attractor through a numer- 
ical experiment in Section 7. 

In the whole article we use the following conventions: If i? is a Banach space, then we 
denote by | • 1^; the norm of E and by (•, •)e{=- \ ■ W) the inner product of E^ if E is 
even a Hilbert space. In the cases E = M" for some n G N or i? = T^ x M^, where 
T^ denotes the two-dimensional torus, we just write | • | = | • |_b and (•, •) = (•, ■)e for 
the norm and inner product, respectively. We also use |2:| to denote the absolute value 
of a scalar z G M or a complex number z G C. Further, by * we denote the complex 
conjugate of a complex number or a complex-valued function. For a separable Hilbert 
space E we denote by C{E) the Banach space of linear, bounded operators form E into 
E equipped with the operator norm | • \c(e) and C2{E) equipped with the inner product 
{T,S)c,^l^E) '■= YlneNi'^^n, Sen) E, T,S £ £2(^)1 denotes the Hilbert space of Hilbert- 
Schmidt operators from E into E, where (en)neN is a orthonormal basis of E. Further, 
we denote by ($7, J^, P) always a probability space and E the expectation w.r.t. P. For 
< p < oo and E a separable Banach space, we write U'(i^,E) and LP{0,) := LP(r2,M) 
if i? = M for the Banach space of i?~valued integrable random variables X : 0, ^ E (in 
fact, equivalence classes of random variables, where X ^ Y ii X = Y P-a.s.) equipped 
with the norm l-^l^^p = E(|X|^). In the same way we define the spaces L^(M) with 
respect to the Lebesgue measure. We denote by 13(E) the Borel cr-algebra of a metric 
space E. Finally, for a multi-index 5 = (5i, . . . , 6n) £ N^ we set \6\ := 6i + ■ ■ ■ + 6n 
and denote the partial derivative operator by D := d' '/dx^^ . . .dx^ and D := id if 
\6\=0. 

2 Motivation for the use of fractional noise 

Since solutions of the Navier-Stokes equations for very turbulent fluids, i.e. at large 
Reynolds numbers, are unstable in view of the sensitive dependence on the initial condi- 
tions that makes the fluid flow irregular both in space and time, a statistical description 
is needed (see e.g. [9]). Based on this, assume that v{x, t), x G M^, t > 0, is a measurable, 
time stationary, space homegeneous and (local) isotropic random fleld on a probability 
space (il, J-", P), such that v satisfies a;-wise, cj G O, the Navier-Stokes equations in two 
dimensions. Due to a phemenological approach, first introduced by Kolmogorov ([18]) 
in three dimensions and by Kraichnan ([19]), Leith ([22]) and Batchelor ([4]) in the two 
dimensional case, we have the following phenomenological correspondence (see e.g. [5]) 
between the relation of the second-order structure function 

E{\v{x + r,t)-v{x,t)\^) =C|r|"-i (2.1) 

and the relation of the energy spectrum E{-) of the velocity field v 

E{k) = Ck-'', (2.2) 



where C,C > are some constants, 1 < a < 3, r G M^ and A: > in the inertial 
subrange. In particular, for a = 5/3 we obtain the famous Kolmogorov's two-thirds law 
and Kolmogorov's five-thirds law (or Kolmogorov energy spectrum), respectively. 
The connection to the fractional Brownian motion gives now Taylor's frozen turbulence 
hypothesis ([44]) which informally assumes that the spatial pattern of turbulent motion 
is unchanged as it is advected by a constant (in space and time) mean velocity V, 

\V\ := (X^j^j ) ! l^t us say along the x axis. Mathematically, Taylor's hypothesis says 
that for any scalar-valued fluid-mechanics variable ^ (e.g. Vi, i = 1,2) we have 

The frozen turbulence hypothesis enables us to express the statistical characteristics of 
the space differences v{x + r,t) — v{x,t) in terms of the time differences v{x,t) — v{x,t + s) 
corresponding to a fixed time t. Indeed, (2.3) implies Vi{x,t + s) = Vi{x — Vs,t) and 
therefore by (2.1) we deduce 

E {\v{x, t) - v{x, t + s)|2) = C|F|"-is"~i (2.4) 

in the inertial subrange along the time axis. For a discussion under which conditions the 
frozen turbulence hypothesis is valid, see [27]. 

Note that due to our derivation the properties (2.2) and (2.4) are closely related to each 
other! 

Now comparing (2.4) with the statistical property E [IPl^g ~ /^t P) — ^ of the frac- 
tional Brownian motion (fBm) (/3/^)jg]R with Hurst parameter H G (0, 1) indicates that 
it is reasonable to model the random velocity field v with noise driven by a fBm. With 
this motivation Shao [38], Sreenivasan et al. [17] and Papanicolaou et al. [29] proposed 
models of the velocity field driven by a finite-dimensional fBm, always with special in- 
terest in the case H = 1/3 in view of the Kolmogorov energy spectrum. 
To capture both statistical features, (2.2) and (2.4), we introduce in Section 4 a two- 
dimensional, incompressible, stationary, homogeneous and isotropic random velocity 
field. These assumptions (Gaussian statistics, stationarity and isotropy) are quite com- 
mon in random mathematical models of turbulent fluids, see e.g. [3]. 

3 Stationary fractional Ornstein-Uhlenbeck process 

As will be seen in the next section, the stationary solution of our stochastic evolution 
equation will be given by an infinite series of stationary fractional Ornstein-Uhlenbeck 
processes. For that purpose we recall in this section from [8] basic properties of the 
stationary fractional Ornstein-Uhlenbeck process. 

First, we remind that a real-valued and normalized fractional Brownian motion (fBm) 
on M with Hurst parameter H E (0, 1) is a Gaussian process (3^ = (/3^)tgi8 on a prob- 
ability space (fi, J",P), having the properties /3^ = P-a.s., E (/3/^) = 0, E (/3f /3f ) = 
2 (N^ + \^\ ~ 1^ ~ ■^\ } ) s,t G M, with continuous sample paths P-a.s.. The reader, 



interested in stochastic calculus w.r.t. fBni, is referred to [26]. 
We consider now the fractional Langevin equation 

t 
Xt = Xo-va Xsds + u^VX^I^, t > 0, Xq G M, (3.1) 



where q,A,i^ > 0. The unique stationary solution of (3.1) is given by the fractional 
Ornstein-Uhlenbeck process Yt := u^V^ J_ e~*^*~"^'^"d/3^, t G M, see [8]. The unique- 
ness has to be understood as uniqueness in law in the class of stationary solutions adapted 
to the natural filtration generated by the two-sided fBm /3^. 

The following proposition will be used several times in this article. Here and in the 
following r(x) = J„°° s^~^e~'^ds, x > 0, denotes the gamma function. 

Proposition 3.1. Let Yt, t ^ M, be the unique stationary solution to (3.1) and set 
C{H) := T{2H + 1) sm{T:H)/T^ > 0. 

(i) Then for all t,s £M we have 

oo 

Cov(Yt,Ys) :=E{YtYs) = C{H)^jy / cos{{t - s)iyax)- ^dx. (3.2) 


In particular, for allteR we have Var{Yt) := E (y^^) = \T{2H)H/a^^. 



(a) Fix T > 0. Then there is a constant Ci{H,X,i',a,T) > such that for any 
t, s £ [— T, T] we have 

Ci{H,X,u,a,T)\t-s\''" <E{\Yt-Ys\^). 

(Hi) Then for any 7 G (0,^] there is a constant C2(if, z^, 7) > such that for any 
t, s G M we have 

E{\Yt-Ys\^)<C2{H,u,-f)Xa^''-^''\t-s\^\ (3.3) 

Proof, (i): By Remark 2.4 in [8] we have 

00 

Cov{Yt,Ys) = iy^^X— -^ ^ / e*^* '>'' dx, 

2-K J [va)'^ + x^ 

—00 

where t, s G M and i denotes the imaginary unit. (3.2) follows now by the change of 
variables y = xjiya) and by the fact that sin(-) and cos(-) are odd and even functions, 
respectively. The relation for the variance of Y is the consequence of (3.2) with the 
relations /~ ^dx = ^^i^^™, see (3.241.2) in [15], r(l - F)r(F) = ^^^ and 



T{2H + 1) = T{2H)2H. 

(ii)+(iii): Set C{H) := — ^'^'"^^ and fix without loss of generality — oo < s < t. 

By (i) and the change of variables z = (t — s)vax we have 



oo 

„1-2H 



E {\Yt - y,p) = 2C{H)^^ / (1 - cos((t - s)vax))-—^dx 

(3-4) 

= 2C{H)X{u{t - 8)^" / (1 - cos(z))^- , .. , dz. 

J [(t — s)uay + z^ 


Lower bound: For any —T < s < t <T, where T > is fixed, we obtain by (3.4) 

oo 

E{\Yt-Ys\^)>2C{H)X{u{t-s)f'' inf / (1 - cos(z))-- , ,, , , dz 

^ ^ -T<s<t<T J {{t-s)vaY + z^ 



oo 



2C{H)X{u{t - s)f" / (1 - cos(z))— -^— ^ ^dz. 



,1-2H 

(2rz/a)2 + z2' 



Upper bound: By (3.4) we have for all — oo < s < t 

A-2H 



E {\Yt - y,p) < 2C{H)X{v{t - s)f^ sup /"(l - cos(z)) 

— oo<s<t J 



dz 



((t- s)vaY + z 





= 2C{H)X{u{t - s)Y'' / (1 - cos(z))z-^-"^dz. 



Notice that /g°°(l — cos{z))z^^^'^^ dz = 2 j^ s\v?{z/2)z^^^'^^dz and this indefinite 
integral is finite. Indeed, this follows from the fact that we have {) < H < 1 and 
sm^{z) = I sin(z) - sin(0)p < C{e)z'^'' for any e G [0, 1] and a constant C(e) > 0. 
It is left to prove the relation (3.3) for 7 G (0, H). By (3.2) and using again the Holder 
continuity of sin(-), we get 

00 

E {\Yt - y,|2) = 2C{H)^^ / (1 - cos((t - s)vax))-—^dx 



00 

A /■ x^~2/^ 
= AC{H)^^ sm'^{{t-s)uax/2)- ^dx (3.5) 



00 

f l+2-y~2H 

< C(F,7)Aa2T-2^(zy(t - s)f^ / ^^^ ^dx < 00, 

J 1 + x^ 



for any 7 G (0, i/) with a constant C{H, 7) > 0. D 



4 The model 

The concrete model we study is the following system (M) of equations in non- 
dimensional form 

Tx{t) = v{x{t),t) -x{t), (x(0),i;(0)) gT^ xM^, (Ml) 

v{x, t) = V^^ix, t) = (^(x, t), -^(x, t)) , (M2) 

dipt = uAiptdt + u^Q^dBf, ipo £ V, t>0, (M3) 

where we assume that 
Assumption 4.1. (i) r, z^ > 0. 

(ii) V := {f £ L^'P^''(T^)| fj2 f{x)dx = 0} is the separable Hilbert space with inner 
product {f,g)v '■= ff2 fix)g*{x)dx, f,g £ V, and with orthonormal basis (ONB) 

(6fe("))fcgi^ '■= (^^^'''Ofceii-' where k e K := 27rZ^\{(0, 0)} and i denotes here and 
in the following the imaginary unit. 

(iii) A : ^^{A) C ^ — )• ^ is a linear operator such that there is a strictly positive 
sequence (afc);.^^ C [c,oo) with c > 0, a^ = a-k, Aeu = — afcCfc and a^ — >■ oo for 

|A:| — )• oo. 

(iv) Qa : y — ). y is a bounded linear operator such that there is a positive sequence 
(\/^)fcei^ ^ [0' oo) with yJJ^ = y^Xl^ and Q^Ck = y/Xke-k- 

(v) iyB^\ is an infinite-dimensional fractional Brownian motion in V with Hurst 
parameter H € (0, 1) defined on a probability space ((7, J^, P) by the formal series 

B''(t) = Y,Pk{tyk. tGM, (4.1) 

ki^K 

where [{Pj^ {t))t^^, k £ K^ is a sequence of complex-valued and normalized frac- 
tional Brownian motions, each with the same fixed Hurst parameter H G (0, 1), 
i.e. (3^ = -y=Re{/3^) + iA=I'm{Pj^), where Re{/3^) and Im{/3^) are independent 

real- valued and normalized fractional Brownian motions on M, and different (3^ 
are independent except /3;^^ = {(3^)*- 

So by Assumption 4.1 ^ is a strictly negative, self-adjoint operator with compact 
resolvent and Q^ a positive, self-adjoint operator, respectively. We will show in this 



section that the conditions on the sequences {f3J^)keK, {Xk)keK and {ak)keK in Assump- 
tion 4.1, together with some growth conditions on {Xk)k<aK and {ak)k£K, imply that 
ip{x,t) and the components of v{x,t) = V-^tlj{x,t) are real-valued. 

Under the Assumption 4.1 on the operator A, it follows that vA generates an analytic 
semigroup on V which in the following will be denoted by {St)t>o- Since the spectrum 
of vA is strictly negative and has a lower bound strictly less than zero, the semigroup 



{St)t>o is exponentially stable. In particular, we have |'S'(t)|£(y) < e '^^^^ffes-K^fc < g *c;/ 
for afl t > 0. 

The domain ^(A) of ^ is given by i:'(^) := {/ G V\ EfeeA' «fcK/' efc)yp < oo}. Further, 
we define fractional powers {—A)'^ : D((— A)"^) C F — )• F, 7 > 0, of the strictly positive 
operator {-A) by P((-A)^) := {/ G y| EfceA' "^K/' efc>yP < 00}. P((-^)^) endowed 
with the inner product 

{{-Ayf, {-Ayg)y = ^ al'if, ek)v{g, euYv =■ if, 5)(-A). 

k&K 

for f,g £ V{{—A)'^), becomes a Hilbert space. Notice also that V{uA) = ^^{A) and 

V{{-vAy) = v({-Af) for any 1/ > and 7 > 0. We similarly define V{{-A)^) for 

7 < as the completion of V for the norm | • | (_a)7 • 

We are mainly interested in the special case when A = A, where A denotes the Laplace 

operator with periodic boundary conditions. Then a^ = \k\'^ , k £ K, and D((— A)) = 

VF^'2(T^) n V, where Ty^'^(T^) n V denotes the Sobolev space of functions on T^ whose 

weak derivatives up to order 2 are in V. In particular, we have D((— A)''') = T^^''''^(T^)n 

F for 7 > 0. 

We will only use the following concept of solutions to equation (M3). For that fix T > 0. 

A 13{[0,T]) (g) J^-measurable F-valued process (v^(t))te[o,T] is said to be a mild solution 

of (M3), if for alltG [0,r] 



V'(i) = S{t)ij{0) + S{t- s)Q2dB^{s) (4.2) 



P-a.s., where the stochastic integral on the right hand side of (4.2) is defined by 

t t 

I S{t - s)QhB''{s) := Y, ^ki^"" /"e-(*-^)^"'=d/3f (s) e^, (4.3) 

^=6^ 

provided the infinite series in (4.3) converges in L^(il,y). 

We are mainly interested in strictly stationary solutions of (M3) which are ergodic. 
Therefore, we call a mild solution {ip{t)t>o) strictly stationary, if for all A; G N and for all 
arbitrary positive numbers ti,t2, ■ ■ ■ ,tk, the probability distribution of the F^-valued 
random variable [ip{ti + r), ^^(^2 + r), . . . , ip{tk + r)) does not depend on r > 0, i.e. 

Law {^p{tl + r), ip{t2 + r), . . . , ^'(tfc + r)) = Law (V'(ti), ^(^2), • • • , "0(4)) 

for all ti,t2, ■ ■ . , tk, r > 0. Here Law{-) denotes the probability distribution. We say that 
a strictly stationary mild solution {ip{t))t>Q of (MS) is unique if every strictly stationary 
mild solution of (M3) which is adapted to the natural filtration generated by the two- 
sided infinite-dimensional fractional Brownian motion (4.1) has the same distribution 
as (■i/'(t))i>o. Further, we call a strictly stationary solution ergodic if for all measurable 



functionals p : V ^ M such that E (|p('0(O))l) < oo we have P-a.s. 

T 
lim i /p(V^(t))dt = E(KV'(0))). 


We have the following existence and uniqueness result. 

Theorem 4.2. Suppose Assumption 4-1 holds and assume that there is e > such that 

Y, hc^f'"^ < oo. (4.4) 

keK 

Then there exists a unique ergodic mild solution tp to equation (M3) given by 

t 

Ht) = E V^^"" / e-(*-")'^"'=d/3f (n)efc, t e M. (4.5) 



Proof. The existence of a strictly stationary mild solution which is ergodic follows by 
Theorem 3.1 and Theorem 4.6 in [23]. We only remark that Theorem 3.1 and Theorem 
4.6 in [23] are applicable since we have in particular 

\S{t)Q'2\cHv) = \{-i^Ar{-uA)-''S{t)Q^cHv) = \{-i^AyS{t){-uA)-^Q^cHv) 

< \(-uArSit)\civ)\i-^A)-^QMcHv) = \{-i^ArS{t)\civ) ( E 



A. V' 



k&K ^ "^ / 



for any 7 G [max{0, -ff — e}, 00), a constant C > and any t > 0. Here we used (4.4) and 
the well-known interpolation inequality \{—i'A)"'S{t)\jrfy\ < Ct~"' , see e.g. Theorem 
2.6.3 in [35]. 

Now assume that we have two strictly stationary mild solutions ip and ^p to equation 
(M3). Notice that 



l'0(t) - V(t)|y = \S{t)iiPiO) - V(0))ly < e~*'^'°^*e^°'=lV(0) - V(0)lv' ^0 for t 



CXD 



P-a.s. and this implies uniqueness in the sense of our definition. The representation 
(4.5) of the mild solution is just the consequence of the definition of a mild solution and 
stationarity. D 

The next remark will be useful in this and the following sections. 
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Remark 4.3. Let ip be the unique ergodic mild solution to equation (M3) given by (4.5). 
For any i G M and k £ K we set in the following 



=di?e(/3f )(u) + i^y ^^^ J e-^'-->"''dIm{P^){u) 



oo 

t , t 

2" 




U^ / 6"*^*"'"''''°'=' 



-00 



= : '4>k,Reit) +ill'k,Im{t) 

and therefore ip{x,t) = Yl,k&K'^k{t)e-k{x), t £ ^, x £ T^. Further, observe that by 
Assumption 4.1 we have 



{Mt)r = {^/^ku" / e-^'-^>^''dp^{u)) =v^i^^ / e-(*-")^°'=d(/3f)*(n) 



-00 
t 



= TaI;^^ y e-('-^>^-Hp\{u) = ^_fc(t). 
—00 
In particular, for any s, t G M and k,k' £ K we obtain 

jg/T /,w? . XN*N ^ f2E(V'fc,Re(t)V'fc,fle(s)) = 2E('(/;fc,/m(i)'0fc,/m(s)) if A; = fe' 
[Vk[ )[Vk'[S)) ) 1^ iffc^yt' 

and 

E{\Mt)-Ms)f) = 2E{\lPk,Re{i) - ^kM^)f) = 2E{\lPk,Im{t)-A,Imis)f)- 

Therefore, to compute E(^fc(i)(V'fc'(*))*) o^" IE(|'0A;(t) — V'fc('S)| ) we only have to compute 
the associated real part and multiply it by two. 

Theorem 4.4. Suppose Assumption 4-i holds. Further, assume that there is m G No 
and 7 G (0, 1) such that 

^ Afca27-2^|A;|2- < 00 and J^ Afca^2//|^|2n^+2^ ^ ^^ (4_g) 

Then there is a unique ergodic mild solution t/j to equation (M3). Further, for all 5 G Nq 
with \6\ < m there is a version of D ip (again denoted by D ip) such that 

D^i; G C"(t2 xM) 

F-a.s. for any e G (0,min{7, //}). In particular, D t/j is real-valued. 
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Proof. By (4.6) and since Assumption 4.1 holds Theorem 4.2 implies the existence of a 
unique ergodic mild solution ip to equation (M3). Further, for all t G M and x € T^ 
ipixjt) is real-valued, since (tpkit))* = tjj-k{t) and therefore 



(v(x,t))* = ( ^ 4(t)e^<''">) = J2(Mt)ne^^''^>r 

k&K keK 

= Y, V^-fc(t)e*<-''"'> = Yl Mt)e'<'''''^ = V(x,t). 



keK kGK 

Now let m G No and 6 = {5i,52) G Nq with |(^| = (5i + (52 < m. By (4.6) it is clear that 
the stochastic process {D ip{t))teR, defined by the formal Fourier series 

t 

keK J^ 

is a well-defined 'D((— A)^)-valued stochastic process for some ^ G M. But D ip is in 
general a function (and not only a generalized function) if (^ > 0. However, < (^ < 7 is 
already assured by 

E(|i?V(t)l?_^)c) < E <\^n\m?)\k?^ = mH)H E Xkal^-'^'lkl'^ < 00, 

keK keK 

where we used \D^ek{x)\'^ < |/cp™, Remark 4.3, Proposition 3.1(i) and (4.6). Notice also 
that D^'il){x, t) is real-valued by using the same argument which leads us to conclude that 
ijj{x,t) is real-valued. Again, by \D ek{x)\'^ < \k\'^'^, Remark 4.3, Proposition 3.1(iii) 
and (4.6) we obtain for all t, s G M and x G T^ 



E{\D^iP{x,t) - D^i;{x,s)\^) = ^ H\Mt) - Ms)\^)\D^ekix)\' 

keK 



keK 
< C{H, u, e) J2 Xkaf'^^\k\^'^\t - s\^' < 00 

keK 
for any e G {0,nim{j,H}) and some constant C{H,v,e) > 0. Similarly, using (4.6) and 

\D'ek{x) - D'ek{y)\ < \kr\eu{x) - efc(y)| < C{r^)\kr+^\x - yY< 

for any rj G (0, 1) and a constant C{rj) > 0, we get for all t G M and x, y G T^ 

E{\D^ij{x,t) - D^i;{y,t)\^) = ^ E(l4(i)l')l^'^efc(x) - I)%(y)|2 

keK 

< C{e)T{2H)H ^ Xkaf\k\^'^^^'\x - y\^' < 00 

keK 
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for any e G (0,7) and a constant C(e) > 0. Therefore, we obtain for all t, s G M and 

E{\D^i;{x, t) - D^il^iy, s)\^) < C{H, v, 6, e){\t - s\'^' + \x - yj^^) 

for any e G {0,rain{j, H}) and a constant C{H,v,6,e) > 0. As D ip{x,t) is a normal 
real-valued random variable, we have 

E(|Z)V(x, t) - D'i;{y, s)\^^) < C{H, u, 6, e, n) {\t - sp- + \x- y^-) 

< C{H,v,5,e,n){\t - s\^ + \x - y\y 

for all t, s G M, x, y G T^, n G N, e G (0, min{7, //}) and a constant C(i?, z^, 5, e, n) > 0. 
Theorem 3.4 in [12] implies now that there is a version of D ^fj (again denoted by D tp) 
such that P-a.s. 

D^i; G C"(t2 xM) 

for any e G (0,min{7, //}). D 

Remark 4.5. In [28] Nualart and Viens established an analogue regularity assertion as 
in Theorem 4.4 for the mild solution of the fractional stochastic heat equation on the 
circle, but they did not consider the partial derivatives of the mild solution. 

Corollary 4.6. Suppose all assumptions of Theorem 4-4 hold. In particular, assume 
that there is m G No and 7 G (0, 1) such that (4-6) is satisfied and let ip be the unique 
ergodic mild solution to equation (M3). Then there is a version of ip (again denoted by 
ip) such that P-a.s. 

V'GC7(R,C""(t2)). 

Further, for all —00 < Ti < T2 < 00 there is a positive random variable 
K = K{H,u,m,j,Ti,T2) : n^ [0, 00) withE{K'^) < 00 such thatF-a.s. 

IV'(w)lc([Ti,T2],C'"{T2)) < K{^)- 

Proof. Recall from the proof of Theorem 4.4 that for any 6 G Nq, \5\ < m, e G 
(0,min{7, ff}), fc G N with ek > 1 and s,t G M, 2; = {xi,X2),y = (2/1,2/2) S T^ we 
have 

E{\D^i;{x, t) - D^i,{y, s)\^^) < C{H, u, m, 6, e, k) {\t - s^^^ + |^ _ y|^2fc^ 

< C{H, u,m, 6, e, k){\t - sf' + \xi - yif' + \x2 - 2/2^'') 

for some constants C{H,v,m,5,e,k),C{H,v,m,5,e,k) > 0. By Theorem 1.4.1 in [20] 
{Kolmogorov's continuity theorem for random fields) for / G N, e2l > 3, —00 < Ti < 
T2 < 00 and any 

^ e2l lHf~^\ /e2/-3\ , ^^, 
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there is a positive random variable K = K{H, i^, ?7i, (5, e, /3, /, Ti, T2) : $7 — t- [0, 00) with 

E{K^^) < 00 (4.7) 

and a version of ip (again denoted by Tp) such that for all t,s £ [Ti,T2] and x = 
{xi,X2),y = {yi,y2) £ T^ we have P-a.s. 

\D^i:{x,t){co) - D^^{y,s){uj)\ < K{io){\t - sf + |xi - yif + \x2 - ysl^- 

In particular, for t, to G [Ti,r2] and x = {xi,X2),xo = (2;o,i) 2:0,2) £ T^ we have P-a.s. 

\D^^ix,t){uj)\ < \D^i;{x,t){io) - D^^j{xo,to){uj)\ + \D^i^{xo,to){uj)\ 

< K{uj){\t - tof + |xi - xo,if + \X2 - xcsl'^) + \D^ij{to,xo){uj)\ 

and therefore P-a.s. 

sup sup \D^tl^{x,t){cv)\ < Kioj)Ci/3,Ti,T2) + |DV(^o,io)(^)| (4.8) 

te[Ti,T2]xeT2 

for some constant C(/3,Ti,T2) > 0. The assertions of the corollary follow now by (4.7), 
(4.8) and since we have 

E{\D'i;{xo,to)\') = j; AfcZ.2^E I I / e-(*°-")'^"'=d/3f (^)n |D%(xo)P 

k^K \ -co / 

keK k 

for every xq G T^, to £ I^ and 5 G Nq, \6\ < m, where we used (4.6), Proposition 3.1(i) 
and Remark 4.3. D 

We apply Corollary 4.6 to state an existence and uniqueness result for the transport 
equation (Ml). Consider the differential equation (Ml) as a first order system 

S(:S) = /*M..(M.W. *(*))). (:J^|) = (!) € T^ X M^ (49) 

where r > and f^(uj),T : K x T^ x M^ _^ M^ is defined by 

it,x,y) ^ /^H,.(t,x,y) = (i(v^^( Jt)(^) - y)) " ^"^'^^^ 

Here V'(") ')i^) with w G il denotes a realization of the ergodic mild solution of (MS). 
We say that (Ml) has a unique local C"^ -solution P-a.s. for some m G N if for all 
(x, y) G T^ X R^ there is an open interval /(w) C M including and a function 

I'^j'A e c' (/(a;),C"(T2 x M^T2 x M^)) ^ 

which satisfies uniquely the equation (4.9) for all t G /(w) P-a.s.. We say that (Ml) 
has a unique global C^ -solution P-a.s. for some m G N, if (Ml) has a unique local 
C™-solution with I{lo) = R P-a.s.. 
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Corollary 4.7. Suppose Assumption 4-1 holds. Further, assume that there is m £ N, 
m > 2, and 7 G (0, 1) such that 

j; A,a27-2/f |^|2m < ^ ^^^ j;Afca^2/f|fc|2-+2^<oo. (4.11) 

k&K keK 

Then (Ml) has a unique global C"^~^ -solution F-a.s.. 

Proof. By Corollary 4.6 there is a version of the strictly stationary solution ip to equation 
(M3) (again denoted by V) such that V G C(M,C™(t2)) P-a.s.. This imphes that 
U{uj),T e C(M,C""-1(T2 X M^M^)) ^i^gj.g j^^^^^^ ig defined in (4.10). Therefore, see e.g. 
Appendix B in [1], (Ml) has a unique local C™~^-solution P-a.s.. To establish P-a.s. 
global solutions we have to find locally integrable, positive functions a^^, /S^ : M — )• [0, 00) 
which may depend upon the realization uj £ Q, such that P-a.s. 

l/v>H,r(i,2;,y)| < a^(t)|(x,y)| + I3^{t). 

By Lemma 6.4(ii) (see Section 6 below) there is a constant K{uj) > such that P-a.s. 

|^V(^,t)(a;)|' + \J^^^Pix,t){u)f < m)i^)\liij2) < i\t\ + K{u;)f 
for all t G M. Hence, we have for all t G M, x = (xi, X2) G T'^ and y = (2/1,2/2) G I^^ 
l/,/.H,r(*,2;,y)p = yj + yi + \—tp{x,t){uj) - yi\ + | - —'il^{x,t){uj) - y2\ 

\—^(x,t)iu;)f + \—i;ix,t)iu;)\' 

<3\{xi,X2,yi,y2)\^ + 2{\t\+K{uj)f 
P-a.s. and in particular 

\Ui^lr{t,x,y)\ < V3\{x,y)\ + V2{\t\ + K{u;)) . 

D 

5 Matching desired statistical properties of the velocity field 

We suppose that Assumption 4.1 holds and set ^ = A (and thereby q^ = \k\'^, k G K). 
Further, assume that 

^A,|fcp+4^-4^<oo (5.1) 

k&K 

for some 7 > 0. So by Corollary 4.6 there is a unique ergodic mild solution ^ of (M3) 
and there is a version of ^/J (again denoted by ^) such that P-a.s. if) £ C (M, C"^(T^)). 
Since v = V"*"^, we have 

v{x,t) = (""'J^'Sl = y^ ^f ^f \Mtyk{x) = V Vk{t)ek{x), t G M, X G t2, (5.2) 

\v2{x,t)J ^^ \-kJ ^^ 
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where we set Vk{t) := i(_| )ipk{t)- Since {il)k{t))t(m.-, k £ K, are mean zero Gaussian 
processes, independent, except (V'fc)* = i/^-k, i" is a mean zero Gaussian random field. 
The autocovariance function i? : T^ x T^ x M x M — >■ M^^^ of v is given by 



R{x,y,t,s) = {E{vi{x,t)vj{y,s)))^^..^^ 
kik2 kf 



E 



keK 

where we set 



^2 ^2^1 1 Xk\k\-^^5k{t - s; H, v) eu[x - y), 



T{2H + l)s\.n{TTH) 



6k(t — s; H, v) := / cos((t — s)z/ A; z) ^dz 

IT J 1 + Z^ 



and where we apphed Proposition 3.1(i) and Remark 4.3. Therefore, v is stationary and 
homogeneous. For k & K, the energy of the Fourier mode k is defined by 

£{k) = ]-E{vkmvk{t)r) = ]:E {\vk{t)\') = hk\'E (\Mt)\') = hk\'E (\M0)\' 






T(2}i)ii^ ,,^,2_4H 



where we again apphed Proposition 3.1(i) and Remark 4.3. 

To ensure isotropy, we set Afc := Cd^D) k G K, for a suitable positive function ( : 
[0, do) — )• [0,oo), but keep in mind that (5.1) should still be satisfied. Then the energy 
of a Fourier mode k £ K depends only on the length of k, in that, £{k) := £{k) = £{k') 
whenever k = \k\ = \k'\ for k, k' £ K. In such an isotropic random field of the form (5.2) 
it is customary to define the energy spectrum in terms of total energy in all the Fourier 
modes of the same length k = \k\ by E{k) := #{k S K\ \k\ = k} £{k). Clearly, E{k) 
can be approximated by E{n) ~ Ck£{k,) with some constant C > 0. 
In general the energy spectrum E{-) can be devided in three ranges: For small |A:|, where 
the energy is injected, E{-) increases in |fc| algebraically. For large \k\, where the energy 
dissipates, we just set \k and therefore the energy spectrum to zero {ultraviolet cut-off). 
For intermediate \k\, in the so-called inertial subrange, E{-) decays \n.\k\ algebraically. 
The spectrum {Xk)k<aK of Q can be chosen so that the energy spectrum of v matches 
experimentally observed energy spectra of a turbulent fluid flow, e.g. the Kolmogorov 



spectrum: E{\k\) oc |A;|~3 and therefore £{\k\) ex Xk\k\ a 1^1"^) i-e- Afc oc \k\ 

Obviously, the decay of the spectrum {Xk)keK of Q as |A;| — )• oo determines the regularity 
of the velocity field v and by this also the regularity of the transport equation (Ml). 
From the physical point of view by applying Corollary 4.7 it is clear that for any spectrum 
in the inertial subrange we have a unique global solution to (Ml) due to the ultraviolet 
cut-off. From the mathematical point of view it is interesting to ask whether there is a 
(unique) solution to (Ml) if we match the spectrum {Xk)keK for all modes without the 
cut-off. In view of the Kolmogorov spectrum and (4.6) in Corollary 4.6 we have 
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+4_ff-4H+2m+47 _ i . |-2| r i-f +2m+47 



|^|__+4i._4«+,,n+47 ^ |fc|-^|A;|-t 
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and — I + 2r?T, + 47 < is satisfied for m = 1 and 7 < «• So by Corollary 4.6 
there is a version of ^p (again denoted by ^p) such that P-a.s. ip £ C (M.,C^{T'^)) and 
V = \7 tp S C (M, C(T^,M^)). Hence, by the classic Peano existence theorem there is a 
solution to (Ml), but uniqueness may fail. 

Further, recall that our main motivation to use fractional noise was to match the statis- 
tical property 

E {\v{x, t) - vix, s)|2) ~ C |t - sp^ 
for f , s > 0, X G T^ and some constant C > 0. We have the following result. 

Proposition 5.1. Suppose Assumption 4-1 holds and that there is m €z N and e > 
such that 

J^Afc|A:|2-+^<oo. (5.3) 

k&K 

Then there is a unique ergodic mild solution ip to equation (M3) and there is a ver- 
sion of ip (again denoted by ip) such that P-a.s. ip £ C (M,C™(T^)) and v = V-^ip G 
C (M, C™'1(T2, R2)^ ^ Further, there is a constant C{H, u) > such that for any t,s eR 
and 2; G T^ we have 

E{\v{x,t)-v{x,s)\'^) <CiH,u)\t-s\^" (5.4) 

and for a fixed T > there is a constant C{H, v,T) > such that for any t,s G [—T, T] 
and X G T^ we have 

C{H, I/, T)\t - sp-f^ < E {\v{x, t) - v{x, s)|2) . (5.5) 



Proof. The first assertion follows by Corollary 4.6 in view of (5.3). Further, by (5.2) we 
have 

E {\v{x,t) - v{x, s)|2) = Y, Ifcl'E (jMt) - Ms)\^) |efc(x)|2 

fceif 



Y,\k\'E{\Mt)-Ms)f) 



keK 

for any t, s G M and x G T^. (5.4) and (5.5) follow now by (5.3), Proposition 3.1(ii),(iii) 
and Remark 4.3. D 

Remark 5.2. It should be noted that (5.3) in Proposition 5.1 is very restrictive and not 
satisfied for H £ [g, 1) if tti = 1 and if we use the Kolmogorov spectrum. But in view of 
the ultraviolet cut-off in the region where the energy dissipates, (5.3) is fulfilled for any 
energy spectrum in the inertial subrange. 
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6 The model as random dynamical system and existence of the 
random pullback attractor 

First we recall some required definitions from the theory of random dynamical systems. 
For the general theory of random dynamical systems we refer to the excellent monograph 

[!]• 

In the following {X, d) is a complete separable metric space and 2 denotes the set of 
all subsets of X. Further, for B £ 2 we denote by B the closure of -B in X and by 
B^ := B \ X the complement of B in X. For x € X and B,C £ 2^ we define the 
semidistance by 

dist(x, B) := inf d(x, b) and dist(B, C) := sup inf d(b, c). 

We make the convention d{x, 0) = oo, where denotes the empty set. 

A family (6{t))t£M.) of mappings on Q into itself is called a metric dynamical system and 

is defined by (0,7^, P, {9(t))teR) if it satisfies the following four conditions: 

(i) The mapping (w, t) i— )• 6{t)uj is .F (8) ;B(M)-J^ measurable. 

(ii) 6{0) = id^ = identity map in 0. 

(iii) {6{t))t£R satisfies the flow property, i.e. 9{t + s) = 9{t) o 0(s) for all s,t £M., where 
o denotes the composition. 

(iv) (^(f))tgR is a family of measure preserving transformations, i.e. F{9{t)~^(A)) = 
F{A) for allAeJ=' and t G M, where 9{t)-'^{A) := {u £ Q \ 9{t)uj G A]. 

We say that a metric dynamical system is ergodic if for all ^ G J-" such that 6{t)~^{A) = A 

for aU i G M, we have F{A) G {0, 1}. 

A random dynamical system (RDS) on X over a metric dynamical system 

(ri, J-", P, {9{t))t£M.) with time M is a mapping 

(j) -.R X il. X X ^ X, {t,u},x) >-^ (j){t,u},x), 

with the following properties: 

(i) Measurahility: </> is i3(M) (g) .Fig) i3(X)-i3(X) measurable. 

(ii) Cocycle property: The mappings (j){t,uj) := (j)(t,uj, ■) : X ^ X form a cocycle over 
9, i.e. they satisfy i;^(0, oj) = idx for all a; G fi and (j){t + s, w) = (j){t, 9{s)uj) o 0(s, uo) 
for all s,t G M,a; G ri. Here o denotes the composition. 

We call a RDS cj) continuous or C'^-RDS if (i,x) i— t- cl){t,Lo,x) is continuous for every 
Lo £ Q and we say that a RDS cpis a C'^-RDS, where 1 < A; < oo if for each {t,Lo) £M.xQ, 
the mapping x i— )• <J){t,uj,x) is fc-times difFerentiable w.r.t. x £ X and the derivatives 
are continuous w.r.t. {t,x) £Mx X for each cj G fi. 
Because of the non-autonomous noise dependence of a RDS, generalized concepts of 
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absorption, attraction and invariance of (random) sets have to be introduced. For that 

we also recah some facts from the theory of measurable (closed) random sets, sometimes 

also called measurable multif unctions, see [11]. 

A set valued map D : Q. ^ 2^ taking values in closed subsets of X is said to be 

measurable if for each x G X the map a; i— t- dist{x, D{uj)) is J-'-B{[0, oo)) measurable. In 

this case D is called a closed random set (of X). 

A universe of closed random sets T> (of X) is a system of non-empty closed random 

sets of X, such that T) is closed under inclusion, i.e. if D and D' are non-empty closed 

random sets of X, such that D'{uj) C D{uj), for all uj G Q, and D £ D, then D' € T>. 

Now let (j) he a, RDS and V a universe of closed random sets. 

A closed random set B is called (strictly) 4>-forward invariant if 

^{t, Lo)B{oo) C B{e{t)uj) (</>(t, uj)B{uj) = B{9{t)uj)) 

for ah w G O, t > 0. 

A closed random set B £ D is called D -absorbing if for any D £ D, uj £ Q there exists 

a time tni^) > 0, the so-called absorption time, such that for any t > toi^j) 

<t){t,e{-t)uj)D{e{-t)oo) c B{u). 

A closed random set A £ D with compact values is called random D-attractor of the 
RDS (/) if A is strictly (/(-forward invariant and for any oj £ Q we have 

dist ((p{t,e{-t)u)D{9{-t)u), A{uj)) -;> as t -^ oo 

for any D £ D. 

If there exists a random P-attractor, then the attractor is already unique. Indeed, 

suppose we have two attractors Ai £ D, i = 1, 2. It follows for any a; G 17 that 

dist{Ai{uj),A2{uj)) = lim dist ((j){t, e{-t)uj)Ai{e{-t)u)),A2{uj)) = 0. 

Therefore, Ai{uj) C ^42(0;) for any oj £ Q. Similarly, we can find the contrary inclusion. 
So the 2?-attractor is unique. 

Theorem 6.1. Let (p be a continuous RDS and D a universe of closed random sets. 
In addition, we assume the existence of a (j)-forward invariant and T> -absorbing closed 
random set B with compact values. Then the RDS has a unique random V -attractor 
given by 

A{uj) = Pi (l){t,e{-t)oj)B{9{-t)uj). 

Proof. See Proposition 9.3.2 in [1] or Theorem 2.4 in [37]. D 

Remark 6.2. The existence of a random attractor for a RDS (p implies the existence of 
an invariant measure for (/>, see [10] for the definition and more details. In particular, 
all invariant measures are supported on the random attractor. 
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Now we come back to our model. 

Proposition 6.3. Suppose Assumption 4-1 holds and that there is m G No and'j G (0, 1) 
such that 

^ A,a^^-2^|A;|2™ < oo and Y. Xkaf"\kr+'^ < oo. (6.1) 

fceX k&K 

Then the quadruple {i^,J^,F,{9{t))teM.) defines an ergodic metric dynamical system, 
where 

• Q = C(M, C""(T^)) equipped with the compact open topology given by the com- 
plete metric (i(V', V") := Yl'^=i IV' - V'|n/(2'^(1 + \ip - i'\n)) , where \^ - V'|n := 

SUp_„<t<„ IV'(t) - tjj {t)\ cm (J2), 

• T is the associated Borel a -algebra, which is the trace in Q of the product a -algebra 

(^(C""(T2)))«, 

• P is the distribution of the ergodic mild solution of (M3), 

• (^(t))tg]R is the group of shifts, i.e. 9{t)'ip{s) = ip{t + s) for all t.,s ^M. and ip £ i}. 

Proof. By Corollary 4.6 and assumptions (6.1) the ergodic mild solution ^ to equation 

(M3) is realized on Q := C (M, C™(T^)). Endowing Q, with the compact open topology 

makes Q a Polish space, actually a Frechet space. 

The group of shifts (0(i))tgR on Q defined by 9{t)'ilj{-) = 4^{t + •) for t G M satisfies the 

flow property and is measure preserving, since P is the distribution of the ergodic mild 

solution of the equation (M3). 

Obviously, t i— )• 6{t)jp is continuous for all ip G Cl and ip i— )■ 6{t)ip is continuous for all 

t gR. Therefore, by Lemma 2.1 in [11] (t, -0) i-^ ^(i)V' is B{R) (g) J"- J" measurable. 

Hence, (ri, J^, P, {0(t))t(zM.) defines an ergodic metric dynamical system. D 

The next lemma will be used several times in this section. 

Lemma 6.4. Suppose all assumptions of Proposition 6.3 hold. In particular assume that 
there is m £ Nq and 7 G (0, 1) such that (6.1) is satisfied and let {il,T,F, {6{t))t£M.) be 
the ergodic metric dynamical system introduced in Proposition 6.3. Then the following 
assertions are valid. 

(i) For all —00 < Ti < t < T2 < 00 the mappings ip 1— )• \tp{t)\(jm(j2\ and ip 1— )• 
IV'lc([Ti,T2],c™(T2)) o,re F-B{R.) measurable and there is a constant 
C = C{m,H, z/, 7,Ti,T2) > such that 

E [\idnit)\l.m(j2-^j < K (^Krfnlc([Ti,T2],c™(T2))j ^ ^' (6-2) 

where id^ : il — t- (7, tp ^^ tp. 
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(ii) There is a {6{t))teR-invariant set J^ 3 0.q = Q.Q{m,'^,H,v) C O (i.e. 9{t)^'^Q.Q = 
f^o for all t G MJ with P(r2o) = 1 such that for all 5 > Q, ip ^ f^o there is a constant 
C{'il>) = C{ip,m,^, H,i',6) > such that 

mt)\c^^j2)<S\t\+C{i^) (6.3) 

for all t G M. In particular, the mapping 

\0 fortP^Qo 

with T > is well-defined and J-'-B{M) measurable. 



Proof, (i): Since il is endowed with the compact open topology, the mappings ip i— )• 
\'ip{t)\cm(j2^ and Tp I— )• |V'lc([Ti,r2],c™(T2)) for all — oo < Ti < i < T2 < oo are T-B{R.) 
measurable. Further, by Corollary 4.6 for any —00 < Ti < T2 < 00 there is a positive 
random variable is: (Ti,r2) = K{m,-f,H,u,Ti,T2) G L'^{n,I'^,F) QL^{n,T^,¥), where 
J^^ denotes the completion of J- w.r.t. P such that P-a.s. 

\'4'\c([n,T2],C"^{T^)) <K{Ti,T2,ip), (6.4) 

since K may be measurable only with respect to the completed a-algebra J- . So for all 

— cxD < Ti < t < r2 < 00 we have 



(6.5) 



E[\idn{t)\cm(^^2^) <E[\idn\c{[Ti,T2],C'^iT^))) =^ (Nnlc([Ti,T2],c™(T2))) 

<E^{K{n,T2)) <oo, 

where E is related to the extension of P to J- . 
(ii): We have by (6.4) 

sup |6'(r)V'|c([Ti,r2],c-(T2)) = l'0lc([ri,r2+i],c-(T2)) < K{Ti,T2 + l,tp) (6.6) 

re[0,l] 

and 

E( sup \0{r)idn\c{[Ti,T2],c^{T^))) =^{\'>'dn\c{[Ti,T2+i],C"'iT^))) 

re[o,i] (6.7) 

= E^{\idn\c(iT^,T2+ilc^{T^))) < eF{K{Ti,T2 + 1)) < 00. 

Taking (6.6) and (6.7) into account, Proposition 4.1.3 in [1] {the dichotomy of linear 
growth for stationary processes) with the measurable mapping tp 1— )• |V'Ic([o,i],C""(t2)) 
implies that 

,. l^(0V'lc([0,l],C™(T2)) 
lim sup j-j = 

t^-itoo \t\ 
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on a (0(t))jg]K-invariant set J^ 3 Qq G ^ with P(ilo) = 1- Therefore, for any 5 > 0, 
'0 € f^o there is a constant T{5, ip) = T{6, m,ip) > such that 

|V'(t)|c™(T2) = \9{t)lp{0)\cm(^j2) < |6'(t)V'|c([0,l],C'"(T2)) < S\t\ 

for \t\ > T{5,ip). Hence, by (6.4) for any 5 > 0, ip £ il.() we have 

|V(t)|c'"(T2) = \e{t)^mcm^r2) < 6\t\+K{-T{6,i;),T{6,i;),i;) (6.8) 

for all t G M. Therefore, for r > the mapping 

, ^ fX°ooe^lV'('S)lc™(T2)^^ forV'Gfio 
\o for V ^ ^0 

is well-defined and J^-B{M) measurable, since for r > and n G N the mappings 

//-n^"l^(^)lc'™(T2)^^ for V' e r^o 

^ \o for i^^Qo 

are finite, J^-S(M) measurable and the '(/'■wise limits for n — t- oo are finite by (6.8). D 

Proposition 6.5. Suppose Assumption 4-1 holds and that there is m £N, m >2, and 
7 G (0, 1) such that 



^Xkal'-^"\k\'^"' <oo and ^ Xkaf" \k\'^"'-^^^ < oo . (6.9) 

k&K k&K 

Then the function (/) : M x J] x T^ x M^ — ^ ^^ x M^ 



(t, V, (x, y)) ^ ct>{t, V, (x, y)) := (/>(t, V) T j := T * 

defines a C^^^-RDS over the ergodic metric dynamical system {Q,J-,¥,{9{t))t(^^) in- 
troduced in Proposition 6.3, where {^/.l) G T^ x M^ is the unique global C"^~^ -solution 
for T > 0, il^ e^ and (x, y) G T^ x M^ at time t eR to 

with f^^r : M X T^ X M2 ^ M^ defined by 

{t,x,y) ^ f^At,i^,y)) = [i_^^±^l'^^t^ _y^ 

Here we change fi to $7 := rig, i.e. to the (9(t))t£M.-invariant set Qq introduced in Lemma 
6.4 (ii). 
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Proof. First notice that by (6.9) and Corollary 4.7 for all {x,y) G T^ x M^, V G ^ (i-e. 
ip G f^o) there is a unique global C"^~ ^-solution (i^[J) to equation (6.10). 



Since t i— )■ (/>(t,'(/')(^) is continuous for every (ip,{x,y)) G x T^ x M^ and (x,y) 



1-^ 



(^(t,V')(^) is continuous for every (t, V') G M x i7, to prove the measurability of (p by 
Lemma 1.1 in [11] we only need to prove the measurability of ^ i— >• (j){t,tl^)[^) for every 
(i, {x,y)) G M X T^ X M?. But this measurability is obvious: Since Q is equipped with 
the compact open topology, tp i— )• ip{x,t) is measurable for every (x,t) G T^ x M and 
therefore also ijj i— )• (j){t, ijj) (|^) . 

The cocycle property is just a consequence of the uniqueness of the solution to equation 
(6.10). D 

The next theorem ensures the existence of the random P-attractor. 

Theorem 6.6. Suppose all assumptions of Proposition 6.5 hold. Then the C^^^-RDS 
(j) defined in Proposition 6.5 has a unique random V-attractor A, where the universe of 
closed random sets T> is given by 

V := {D \ D is a closed random set ofT x M with rDiip) •= sup |(x,y)| < oo 

(x,y)eD(V) 
ct 



and r£){0[—t)Tp)e '^* — )■ for t — t- oo and any c > 0, tp £ ilj 
Further, for any 5 > 



(6.11) 



?5 



B'm := < 





r2 ., m2| |„,|2 ^ (1 + '^) / „i^L/,/„.M2 



{x,y) G T^ X R'\ \y\' < '—-^ / er\^p{u)\(.^^J2^du 



-oo 



> , ip gQ, (6.12) 



is a T>-ahsorhing and (j)-forward invariant closed random set. In particular, for any 
6 > we have 

^(V) = fl ^{t,e{-t)i;)B'{e{-t)i;), i;en. 
teN 

Proof. In view of Theorem 6.1 we only have to show that for any 6 > the random 

set B^tjj) := B (t/j) defined in (6.12) is a P-absorbing and (/>-forward invariant closed 

random set. 

Notice here that ip i— )• J_^ e^ \ip{u)\'^i ,j2)du is measurable and finite for every ^ G il. 

This follows by Lemma 6.4(ii) and recall here again the change of il in Proposition 

6.5. Further, the random set valued map B takes values in closed bounded subsets of 

T^ X M^ and is measurable since the random radius ip i— )• rc{ip) = svLpu^y\^(ji^\ |(x, y)\ is 

measurable, where C is a bounded closed random set. Therefore, i? is a closed random 

set. 

In the following, for any set C C T^ x M^ we set 

7rT2(C) := {■Kf2{x,y)\{x,y) G C}, vrB;2(C) := {■K^2{x,y)\{x,y) G C}, 
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where 

7rTr2 : T^ X M^ -> T^ {x, y) ^ x, and 7riR,2 : T^ x M^ -> M^ (x, y) ^ y. 

Next we prove that f/^ i— )• -B('0) is (/)-forward invariant, i.e. (l){t,ij))B{'il)) C B{9{t)ip) for 
all t > 0, "0 G ^- So to prove the i;^-forward invariance of B we only have to show that 

sup |7riR,2(0(t,V')( ))|^= sup {|x(t)p I x(0) = i/} < sup |yp 



(x,y)eiJ(V) Vyy (:r,y)eB{V') (x,y)&B{e{t)i>) 

for any t > 0, V G f^, since TTf2{B{9{t)ilj)) = T^ for any t > and ip £ Q. 
To estimate |i;(t)P for t > 0, we take the inner product of the equation 

Tx{s) = V 'iIj{x{s),s) — x{s) 

with x{s) and obtain 

r^^|x(5)|' = {x{s),V^^l^ix{s),s)) - \x{s)\^ < l|V^V(x(s), s)|2 + l|x(s)|2 - \x{s)\^ 

where we used {zi,Z2) < ^l-ZiP + ^|-Z2p for zi,Z2 G M^. 
Multiplying by -e^ on each side gives 

^(e-|x(s)|2) < -e^|V^V(a;(s),s)|2. (6.13) 

as r 

By integrating the inequality (6.13) from to t we get 



so that 



t 
er\x{t)\'^ - |x(0)|2 < ^ f er\V^^p{x{s),s)\'^ds 





t 
|x(t)|2 < e~^\±{0)\^ + - I e-^|V^V(a;(s), s)pds 



t 
< e-r\x{0)f + - j e-^|V(s)|^i(Tr2)(is 





"~= e"^|x(0)|^ + - / er\il}{u + t)\'^(ji(j2^du. 
-t 



Therefore 





\TT^2m,^)l^}^^^\)\^ ^ (6.14) 
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Now by (6.14) and the definition of B(^) we obtain for t > 



sup \TT^2{(j){t,7p)i ])f<e~T sup |yp + - / e^|V'(w + i)lci(T2) 

y)&BM) \yj (x,y)eBm T J 



{x,y)&B(i,) \y/ {x,y)eB{i,) _^ 





<- — '■ — -e~r I er\'il;{u)\'^i(j2)du-\ er\'ip{u + t)\^,i(j2^du 



— oo —t 

(1 + -^) 7 ^1 // , .m2 ^ , (1 + '^) f ^l u , .m2 



e^\ilj{u + t)\^i^j2)du-\ / er\^l!{u + t)\(-,l^J2^du 



T J -- v.- / 7" 

—oo —t 





(1 + O) / Hi,/ n |2 

= — - — / er\il){u + t)\^^(r^2)du 

— oo 

= sup \y\^. 
(x,j/)es(0(t)t/.) 

So "0 I— ^ -B(^) is (/>-forward invariant. 

Finally, we prove that ip i— )• B[%1)) is D-absorbing. For any D £ D where 2? is defined in 

(6.11) we have 

7rr2{cl){t,e{-t)i;)Dm-t)W) Q 7rr2{B{^j))=T' 

and by (6.14) 



sup \TT^2{<P{t,9{-t)^)r))f 

(x,y)eD{e{-t)'ii>) \yj 







* ■ |2 , -*- / -I / / I + +m2 



<e - sup |y| + - I e^\'4}{u + t-t)\(j^,j2-.du , . 

(x,y)<^D{B{-t)^) T J ' iD.iOj 



e -^ sup |yP + - er\'il;{u)\'^i,j2)du 

(x.v)eD(e(-t)^) T J 



{x,y)eD{e{-t)i,) _^ 



for any t > and ^ E 0. By the definition of the set P the first term on the right 
hand side of (6.15) converges for t — )• oo to zero. The second term tends for i — )• c« to 
^ /_^ er \iIj[u)\'^i ,j2)du. Further, Lemma 6.4(ii) implies that 



-te(l + <^) 



T 

— oo 



e^ |0(u)|(^i/--][.2\dti — )■ as i — )• c« 



for any c, 5 > and ^ G il.. This ensures B G T>, i.e. i? is also P-absorbing. The 
assertion follows now by Theorem 6.1. D 
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Remark 6.7. The universe of closed random sets D defined in Theorem 6.6 contains only 
sets which do not grow with exponential speed (or grow sub-exponentially fast) . Notice 
also that limt_>oo e~^^'f£)(6{—t)^p) = for any c > 0, ip £ i^ where D is a closed random 
set is equivalent to limi_j.oo log(max{l, r£){6{—t)'>p)})/t = for any ^ £ il. Such subsets 
are also called tempered. 

The next step would be to derive (upper) bounds of the Hausdorff dimension of 
the random 2?-attractor which is due to ergodicity of the underlying metric dynamical 
system {Q,T,F,9{t))t£M.) IP-a.s. constant. This will be studied also in the fractional 
noise case in a forthcoming paper. 

7 Simulation 

In this section we visualize the long-time behaviour of particle motions. We simulate 
the motion of lO'^ particles uniformly distributed on the unit square at time zero with 
zero velocities for 5000 time units according to the system (M) with H = 1/3, u = 10~^ 
and the Kolmogorov spectrum, but with different values of r. As already described in 
Section 5, we set the spectrum {\k)k&K of Q for \k\ > 2ttR with some fixed R G N to 
zero, actually i? = 2 in our simulation. Then 



v{x,t)= Yl M t )V'fc(t)e^<''">, xGT^tG 

k&K, |fc|<27ri? 



where ^fe(t) = \/Afci/^ /_^ e~*^*~"^^°'=(i/3^(n). To simulate V'fe) we use the Cholesky 
method (see e.g. [2] pp. 311), since we know by Proposition 3.1(i) and Remark 4.3 the 
covariance structure of V'fc- Given a realisation ofipk^kGK, \k\ < 2ttR, the velocity field 
V can be computed efficiently using the fast Fourier transform algorithm, see Chapter 
4 in [39]. Finally, we integrate the trajectory of the inertial particle using the classic 
fourth-order Runge-Kutta scheme with time step double of that used for the velocity 
field. 

Figure 1 shows the final position of the particles in the phase space associated to four 
different Stokes' number in form of r. The clustering is distinctive for r = 1 and r = 
10^^, but there is almost no clustering for very low and large values of r, i.e. for r = 10~^ 
and r = lO'^ in the experiment. Repeating the simulation under same conditions but 
with different Hurst parameter produces similar results. Therefore, we conclude that 
also the generalized inodel with fractional noise captures the clustering phenomenon 
of preferential concentration. We leave the numerical study of the dependence of the 
clustering on the parameters H and v for a future publication. 
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Figure 1: Snapshots of the final positions of 10^ particles in the phase space associated 
to different Stokes' numbers r. a.) r = 10^^, b.) r = lO^"*^, c.) r = 1, d.) 
r = lO^. 
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